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Abstract. Assume that two robots are located at the centre of a unit disk. Their goal is to evacuate 
from the disk through an exit at an unknown location on the boundary of the disk. At any time the 
robots can move anywhere they choose on the disk, independently of each other, with maximum speed 
1. The robots can cooperate by exchanging information whenever they meet. We study algorithms for 
the two robots to minimize the evacuation time : the time when both robots reach the exit. 

In [9] the authors gave an algorithm defining trajectories for the two robots yielding evacuation time 
at most 5.740 and also proved that any algorithm has evacuation time at least 3 + f + V2 5.199. We 
improve both the upper and lower bounds on the evacuation time of a unit disk. Namely, we present 
a new non-trivial algorithm whose evacuation time is at most 5.628 and show that any algorithm has 
evacuation time at least 3 + ^ + \/3 ~ 5.255. To achieve the upper bound, we designed an algorithm 
which non-intuitively proposes a forced meeting between the two robots, even if the exit has not been 
found by either of them. 


1 Introduction 

The goal of traditional search problems is to find an object which is located in a specific domain. 
This subject of research has a long history and there is a plethora of models investigated in the 
mathematical and theoretical computer science literature with emphasis on probabilistic search 
in [16], game theoretic applications in [3], cops and robbers in [8], classical pursuit and evasion 
in [15], search problems and group testing in [1], and many more. 

In this paper, we investigate the problem of searching for a stationary point target called an 
exit at an unknown location using two robots. This type of collaborative search is advantageous in 
that it reduces the required search time by distributing the search effort between the two robots. 
In previous work on collaborative search, the goal has generally been to minimize the time taken 
by the first robot to find the object of the search. In contrast, in this work, we are interested in 
minimizing the time when the last robot finds the exit. In particular, suppose two robots are in the 
interior of a region with a single exit. The robots need to evacuate the region but the location of 
the exit is unknown to them. The robots can cooperate to search for the exit, but it is not enough 
for one robot to find the exit, we require both robots to reach the exit as soon as possible. 

* This work was partially supported by NSERC grants 
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We study the problem of two robots that start at the same time at the centre of a unit disk 
and attempt to reach an exit placed at an unknown location on the boundary of the disk. At any 
time the robots can move anywhere they choose within the disk. Indeed, they can take short-cuts 
by moving in the interior of the disk if desired. We assume their maximum speed is 1. The robots 
can communicate with each other only if they are at the same point at the same time: we call this 
communication model face-to-face communication. Our goal is to schedule the trajectories of the 
robots so as to minimize the evacuation time , which is the time it takes both robots to reach the 
exit (for the worst case location of the exit). 

1.1 Related work 

The most related work to ours is [9], where the evacuation problem for a set of robots all starting 
from the centre of a unit disk was introduced and studied. Two communication models are intro¬ 
duced in [9]. In the wireless model, the two robots can communicate at any time regardless of their 
locations. In particular, a robot that finds the exit can immediately communicate its location to 
the other robot. The other model is called the non-wireless or local model in [9], and is the same as 
our face-to-face model: two robots can only communicate when they are face to face, that is, they 
are at the same point location at the same time. In [9], for the case of 2 robots, an algorithm with 
evacuation time 1 + yy + a/3 ~ 4.826 is given for the wireless model; this is shown to be optimal. 
For the face-to-face model, they prove an upper bound of 5.740 and a lower bound of 5.199 on the 
evacuation time. 

Baeza-Yates et al posed the question of minimizing the worst-case trajectory of a single robot 
searching for a target point at an unknown location in the plane [4]. This was generalized to 
multiple robots in [14], and more recently has been studied in [11,13]. However, in these papers, the 
robots cannot communicate, and moreover, the objective is for the first robot to find the target. 
Two seminal and influential papers (that appeared almost at the same time) on probabilistic search 
are [5], and [6] and concern minimizing the expected time for the robot to find the target. Useful 
surveys on search theory can also be found in [7] and [10]. In addition, the latter citation has an 
interesting classification of search problems by search objectives, distribution of effort, point target 
(stationary, large, moving), two-sided search, etc. The evacuation problem considered in our paper 
is related to searching on a line, in that we are searching on the boundary of a disk but with the 
additional ability to make short-cuts in order to enable the robots to meet sooner and thus evacuate 
faster. 

Our problem is also related to the rendezvous problem and the problem of gathering [2,12]. 
Indeed our problem can be seen as a version of a rendezvous problem for three robots, where one 
of them remains stationary. 

1.2 Preliminaries and notation 

We assume that two robots R\ and i?2 are initially at the center of a disk with radius 1, and that 
there is an exit at some location X on the boundary of the disk. The robots do not know X, but 
do know each other’s algorithms. The robots move at a speed subject to a maximum speed, say 
1. They cannot communicate except if they are at the same location at the same time. Finally, 
both robots are equipped with deterministic processors that can numerically solve trigonometric 
equations, and as such they are assumed to have the required memory. The evacuation problem is 
to define trajectories for the two robots that minimize the evacuation time. 


For two points A and B on the unit circle, the length of an arc AB is denoted by AB, while 
the length of the corresponding chord (line segment) will be denoted by AB (arcs on the circle 
are always read clockwise, i.e., arc AB together with arc BA cover the whole circle). By ABC we 
denote the angle at B in the triangle ABC. Finally by Aft we denote the vector with tail A and 
tip B. 

1.3 Outline and results of the paper 

In [9] an algorithm is given defining a trajectory for two robots in the face-to-face communication 
model with evacuation time 5.740 and it is also proved that any such algorithm has evacuation 
time at least 3 + | + \[2 > 5.199. 

Our main contribution in this paper is to improve both the upper and lower bounds on the 
evacuation time. Namely, we give a new algorithm whose evacuation time is at most 5.628 (see 
Section 2) and also prove that any algorithm has evacuation time at least 3 + | + y/3 > 5.255 
(see Section 3). To prove our lower bound on the disk, we first give tight bounds for the problem 
of evacuating a regular hexagon where the exit is placed at an unknown vertex. We observe that, 
surprisingly, in our optimal evacuation algorithm for the hexagon, the two robots are forced to meet 
after visiting a subset of vertices, even if an exit has not been found at that time. We use the idea 
of such a forced meeting in the design of our disk evacuation algorithm. 

2 Evacuation Algorithms 

In this section we give two new evacuation algorithms for two robots in the face-to-face model that 
take evacuation time approximately 5.644 and 5.628 respectively. We begin by presenting Algorithm 
A proposed by [9] which has been shown to have evacuation time 5.740. Our goal is to understand 
the worst possible configuration for this algorithm, and subsequently to modify it accordingly so 
as to improve its performance. 

All the algorithms we present follow the same general structure: The two robots R\ and R 2 start 
by moving together to an arbitrary point A on the boundary of the disk. Subsequently R\ explores 
the arc A'A, where A' is the antipodal point of A , by moving along some trajectory defined by the 
algorithm. At the same time, i ?2 explores the arc AA', following a trajectory that is the reflection 
of R\s trajectory. If either of the robots finds the exit, it immediately uses the Meeting Protocol 
defined below to meet the other robot (note that the other robot has not yet found the exit and 
hence keeps exploring). After meeting, the two robots travel together on the shortest path to the 
exit, thereby completing the evacuation. At all times, the two robots travel at unit speed. Without 
loss of generality, we assume that R\ finds the exit and then catches R 2 for our analysis. 

Meeting Protocol for R\ : If at any time to R\ finds the exit at point X, it computes the shortest 
additional time t such that R 2 , after traveling distance to + t, is located at point M satisfying 
XM = t. Robot Ri moves along the segment XM . At time to + t the two robots meet at M and 
traverse directly back to the exit at X incurring total time cost to + 2t. 

2.1 Evacuation Algorithm A, of [9] 

We proceed by describing the trajectories of the two robots in Algorithm A. As mentioned above, 
both robots start from the centre O of the disk and move together to an arbitrary position A on 



the boundary of the disk. i ?2 then moves clockwise along the boundary of the disk up to distance 
7 r, see left-hand side of Figure 1, and robot R\ moves counter clockwise on the trajectory which is 
a reflection of R^s trajectory with respect to the line passing through O and A. When R\ finds the 
exit, it invokes the meeting protocol in order to meet i? 2 , after which the evacuation is completed. 




Fig. 1 . Evacuation Algorithm A with exit position E. The trajectory of robot R 2 is depicted on the left. The movement 
paths of robots R±, R 2 are shown on the right, till the moment they meet at point B on the circle. 


The meeting-protocol trajectory of R\ in Algorithm A is depicted in the right-hand side of 
Figure 1. Clearly, for the two robots to meet, we must have AB = EA + EB. Next we want to 
analyze the performance of the algorithm, with respect to a; := AB, i.e. the length x of the arc that 
i ?2 travels, before it_is met by R\. We also set f(x) := EB. 

It follows that EA = x — f(x), and since 2 sin (EB/2) = EB we conclude that 


/(*) 
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where z is a solution of the equation z — 2 sin 



(i) 


In other words, f(x) is the length of interval EB that R\ needs to travel in the interior of the disk 
after locating the exit at E, to meet R 2 at point B. 

Then, the cost of Algorithm A , given that the two robots meet at time x after they together 
reached the boundary of the disk at A, is 1 + x + f(x). Given that distance x — f(x) traveled by R\ 
until finding the exit is between 0 and 7 r, it directly follows that x can take any value between 0 and 
7r as well. Hence, the worst case performance of Algorithm A is determined by sup xG [ 0 ^{x + f(x)} 
The next lemma, along with its proof, follows from [9]. 

Lemma 1. Expression F{x) := x + f(x) attains its supremum at xq ~ 2.85344 (which is ~ 
0.90827977^). In particular, F(x) is strictly increasing when x G [0,£o] an d strictly decreasing when 
x G [xq, 7r] 


Proof (Sketch). The behavior of F(x), as x ranges in [0,7r], is shown in Figure 2. 


By Lemma 1 , the evacuation time of Algorithm A is 1 + xq + f(x 0 ) < 5.740. The worst case is 
attained for xq — f(x 0 ) ~ 0.30877. 


2.2 New evacuation algorithm 2?(%, <f>) 

We now show how to improve the previously described algorithm and obtain evacuation time at 
most 5.644. The main idea for improving Algorithm A is to change the trajectory of the robots 
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Fig. 2. The performance of Algorithm A as a function of the meeting points of the robots. 


when the distance traveled on the boundary of the disk approaches the critical value xo of Lemma 1. 
Informally, robot i?2 could meet R\ earlier if it makes a linear detour inside the interior of the disk 
towards R\ a little before traversing distance x$. 

We describe a generic family of algorithms that realizes this idea. The specific trajectory of each 
algorithm is determined by two parameters x an d 0 where x £ [tt/2,xo] and </> E [0, f{x)/ 2], whose 
optimal values will be determined later. For ease of exposition, we assume R\ finds the exit. The 
trajectory of i?2 (assuming it has not yet met R\) is partitioned into four phases that we call the 
deployment , pre-detour , detour and post-detour phases. The description of the phases rely on the 
left-hand side of Figure 3. 

Algorithm £>(x,</>)(with a linear detour). R 2 S trajectory until it meets R\ is described below: 

* Deployment phase: Robot R 2 starts from the centre O of the disk and moves to an arbitrary 
position A on the boundary of the disk. 

* Pre-detour phase: R 2 moves clockwise along the boundary of the disk until having explored an 
arc of length y. 

* Detour phase: Let D be the reflection of B with respect to AA f (where A' is the antipodal point 
of A). Then, R 2 moves on a straight line towards the interior of the disk and towards the side where 
O lies, forming an angle of <j) with line BD , until R 2 touches line AA' at point C. From C it follows 
a straight line segment to B. Note that C is indeed in the interior of the line segment A A' by the 
restrictions on <j). 

* Post-detour phase: Robot R 2 continues moving clockwise on the arc BA'. 

At the same time R\ follows a trajectory that is the reflection of R^s trajectory along the line 
AA'. When at time to it finds the exit, it follows the Meeting Protocol defined above. 

Notably, the two robots may meet at point C without having located the exit. Next we consider 
three cases as to where R 2 can be caught by R\ while moving on its trajectory (after R\ has located 
the exit). For all three cases, the reader can consult the right-hand side of Figure 3. As the time 
needed for the deployment phase is independent of where the exit is located, we ignore this extra 
cost of 1 during the case distinction. 


Case 1: i ?2 is caught during its pre-detour phase: The meeting point is anywhere on the arc AB. 
Recall that y < xo, so by Lemma 1 the location F of the exit on the arc FA that maximizes 
the cost of £>(x, </>) is the one at at distance x — fix) from A (see right-hand side of Figure 3). 
The cost then is AB + BF = FA + 2 BF — x + fix)- 












Fig. 3. Illustrations for Algorithm £>(y, 0). 


Case 2: i ?2 is caught during its detour phase: Let G be the point on BC where the robots meet. 
Further, let E be the position of the exit on the arc A'A, and let y := EA. In the following, 
h(y) := EG denotes the length of the trajectory of R\ in its attempt to catch i ?2 after it finds 
the exit. Also, q(y) := BG denotes the distance that travels on BC till it is caught by R\. 
Note that the functions h and q also depend on x and 0; however, while those are fixed, y varies 
with the position of the exit. Lemma 2 below states that h(y ) and q(y) are well defined. 

Lemma 2. y E [x ~ /(x)>x] if and only if the meeting point G of the robots is on the line 
segment BC. Moreover, robot i ?2 can be caught by R\ only while moving from B to C. 

Proof. If the exit is located at point F (i.e. y — X~ /(x))> then the meeting point is B. We 
also observe that if the exit-position coincides with D (i.e. y = %), then the meeting point is 
C. Recall that < f{x)/ 2 = DF/2, and hence point C is in the interior or on the boundary 
of the triangle FDB as it is depicted in right-hand side of Figure 3. Therefore, after time x> 
robot i?2 would approach the exit if it was located anywhere on the arc DE. In particular, if 
y E [x~ /(x)? x] then the meeting point for Algorithm A would be on the arc BA'. In Algorithm 
£>(x,</>), R 2 has a trajectory that brings it closer to the exit. This guarantees that the meeting 
point G always exists, and it lies in the line segment BC (we will soon derive a closed formula 
relating EA and BG). The previous argument guarantees that h(y) is continuous and strictly 
decreasing in y. Notice that EA + EG — AB + BG (since the two robots start from the same 
position A ), which means that 

y + h(y) = x + q(y) (2) 

Hence, q(y) — y + h(y) — y, and as already explained, q(x ~ /(x)) — 0 and q(x) — BC. By 
the mean value theorem, all values between 0 and BC are attainable for q(y) and are attained 
while y ranges in [x ~ /(x)> x]- ■ 

We conclude that if the exit is located at point E , then the cost of the algorithm is y + 2h(y). 
Hence, in case 2, the cost of the algorithm is at most sup y e[ x -f(x),x\{y 2^(y)}- We emphasize 

again that h(y) and q(y) also depend on the fixed parameters x an d 0- 
Case 3: i?2 is caught during its post-detour phase: Clearly, in this case the exit lies in the interior 
of the arc A'D or coincides with A'. At time td — x + 2g(y), robots R\ and i ?2 are located 
at points D and F>, respectively. Then they move towards each other on the arc BD till R\ 
finds the exit. Note that, since DB/2 — sin(y), we have q(x) — sin (y) / cos ((f)). Clearly, the 
closer the exit is to D is, the higher is the cost of the evacuation algorithm. In the limit (as 









the position of the exit approaches D ), the cost of case 3 approaches td plus the time it takes 
Ri to catch i ?2 if the exit was located at D , and if they started moving from points D and B 
respectively. Let G' be the meeting point on the arc BD in this case, i.e. DG' = BG’. We define 
p{x ) to be the distance that R\ needs to travel in the interior of the disk to catch i? 2 , if the exit 
is located at distance x from A. Clearly 6 


p(x ) 


unique z satisfying z = 2 sin 



(3) 


Note also that DG' = p(x) so that the total cost in this case is at most 

U + Mx) = X + 2 sin (x) / COS ((f)) + 2 p(x). 


The following two lemmata summarize the above analysis and express h(y) in explicit form (in 
dependence of x and 0)> respectively. 

Lemma 3. The evacuation time of Algorithm 13 (x, <f>) is 


l + ma x<x + /(x), sup {y+2h(y)}, X + 2 s in (x) / cos (</>) + 2p(x) > , (4) 

l y£[x~f(x),x] J 

where h(y ) (that also depends on the choice of Xi4>) denotes the time that a robot needs from the 
moment it finds the exit till it meets the other robot when following the meeting protocol. 


Lemma 4. For every x > 0 and for every x ~ fix) < V < X; the distance h(y ) that R\ travels 
from A until finding R 2 when following the meeting protocol in Algorithm B(x^4 > ) is 


h(y) = 


2 + (x - y) 2 - 2 cos(x + y) + 2(x - y) (sin(^> + y) - sin (cf> - x)) 


2 (x-y- sin ((f) - x) + sinO + y)) 
In particular, h(y ) is strictly decreasing for 0 < < f{x)/ 2. 


Proof. We start by making some handy observations. For this we rely on Figure 4 (that is a 
continuation of Figure 3). Let H be the point that is symmetric to E with respect to ^L4'. Denote 
with L the projection of H onto the supporting line of DB. Set 6 := BHL , and observe the following 
equation for 6. 
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2 2 2 
7T _ 2 (?r - x) + x - y 
2 2 


7 r BD + DE 
2 2 
X±y_ _ 7T 
2 2 


(5) 


Our goal is to compute h(y) = EG. For this we see 


that E& = eA + lAs + B(5, and therefore 


EG 2 = EH 2 + HB 2 + BG 2 + 2 (e£-H^ + E^-B^ + H^- B&) . (6) 


6 Uniqueness of the root of the equation defining p(x) is an easy exercise. 














Fig. 4. The analysis of Algorithm B(x, </>)• 


2 sin (y) , HB = 2 sin ) > ^G = g(y), and 

eS-H^ = EHHB cos (s/tb) 
k3 -B($ = EH BG cos (g&l) 
H%-B& = TTBBG cos (tt - G-Bif) . 


cos (^EHB^J = cos ( 7 r /2 + 0) = cos 
cos (GBpJ = cos ( 7 t — <fi) = cos (0) 

COS ^7T — GBH^j = COS ^7T + ^ + 0 — = — cos ^ ^ — 4 > 

Substituting the above in (6), we obtain an equation between h{y), q(y ), y, x> and </>• In the latter 
equation we can substitute q(y) using the meeting condition (2), obtaining this way the required 
closed formula for h(y). ■ 

The first natural attempt in order to beat Algorithm A would be to consider £>(x, 0 ), i.e. make 
BC perpendicular to A A' in Figure 3. In light of Lemma 4 and using Lemma 3, we state the 
following claim to build some intuition for our next, improved, algorithm. 

Claim 1 The performance of algorithm £>(x, 0) is optimized when X — Xo ~ 2.62359, and its 
cost is 1 + 4.644 = 5.644. The location of the exit inducing the worst case for £>(x, 0) is when 
y — EA ~ 0.8377T. The meeting point of the two robots takes place at point G (see Figure 3, and set 
<f> — 0), where q(y ) = BG ~ 0.117 ~ 0.236 BC. In particular, the cost of the algorithm, if the meeting 
point of the robots is during i?2 ’s pre-detour, detour and post-detour phase, is (approximately) 5.621, 
5.644 an d 5.644 respectively. 

Note that xo of Claim 1 is strictly smaller than xo of Lemma 1. In other words, the previous 
claim is in coordination with our intuition that if the robots moved towards the interior of the disk 
a little before the critical position of the meeting point xq of Algorithm A , then the cost of the 
algorithm could be improved. 
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2.3 New evacuation algorithm C(%, 0, A) 


Claim 1 is instructive for the following reason. Note that the worst meeting point G for Algorithm 
B(x o 5 0) satisfies BG « 0.236 BC. This suggests that if we consider algorithm £>(xo 5 0) instead, 
where </> > 0, then we would be able to improve the cost if the meeting point happened during the 
detour phase of i? 2 - On one hand, this further suggests that we can decrease the detour position 
Xo (note that the increasing in x cost X + /(x) is always a lower bound to the performance of our 
algorithms when x < #o)- On the other hand, that would have a greater impact on the cost when 
the meeting point is in the post-detour phase of i? 2 , as in this case the cost of moving from B to C 
and back to B would be 2 sin (x) / cos (</>) instead of just 2 sin (y). A compromise to this would be 
to follow the linear detour trajectory of R 2 in B(xo^) only U P to a certain threshold-distance A, 
after which the robot should reach the diameter segment AA! along a linear segment perpendicular 
to segment AA! then return to the detour point B along a linear segment. Thus the detour forms a 
triangle. This in fact completes the high level description of Algorithm C(x, </>, A) that we formally 
describe below. 

In that direction, we fix y,0, A, with x G [7r/2,xo],</> G [0,/(x)/2] and A G [0, sin(x) / cos (</>)]. 
As before, we only describe the trajectory of robot i?2- The meeting protocol that R\ follows once 
it finds the exit is the same as for Algorithms A and £>(x, (/))• 

The trajectory of robot R 2 (that has neither found the exit nor met R\ yet) can be partitioned 
into roughly the same four phases as for Algorithm £>(x,</>); so we again call them deployment , 
pre-detour , detour and post-detour phases. The description of the phases refers to the left-hand 
side of Figure 5, which is a partial modification of Figure 3. 

Algorithm C(x,0, A) (with a triangular detour). The phases of robot R^s trajectory are: 

* Deployment phase: Same as in Algorithm £>(x,</>)- At time 1, R 2 is at point A. 

* Pre-detour phase: Same as in Algorithm £>(x, </>)• In additional time x> R 2 is in point B. 

* Detour phase: This phase is further split into three subphases. 

o Subphase-1: Up to additional time A, R 2 moves along a line segment exactly as in the detour 
phase of Algorithm £>(x, 4>)- Let G be the position of the robot at the end of this phase. 

o Subphase-2: Let C be the projection of G onto AA!. R 2 follows line segment GC till it reaches 
point C. 

o Subphase-3 (Recovering phase): Robot follows line segment CB back to^point B. 

* Post-detour phase: Same as in Algorithm £>(x,</>)- After additional time EA’, R 2 reaches point 
A'. 

At the same time R\ follows a trajectory that is the reflection of R 2 S trajectory along the line 
AA!. If a robot finds the exit, it follows the meeting protocol defined earlier. 

Obviously, Algorithm C (x? 0? cos(0) ) identical to Algorithm £>(x,</>). Moreover, as before 
robots may meet at point C without having located the exit. 

Notice that an immediate consequence of the definition of C(x,</>, A) is that if robot R\ finds 
the exit and meets R 2 during its detour subphase-2 in some point K (as in the right-hand side of 
Figure 5), then 

EA + EA = AB + SG + GK. (7) 

When R\ finds the exit somewhere on the arc A f A, it catches R 2 on its trajectory so that they 
return together to the exit. Note that since robots meet at point C, if the exit is not in the arc 
DB , it is impossible for a robot to be caught by the other robot in subphase-3 of its detour phase. 





c 



Fig. 5. Illustrations for Algorithm C(y, 0, A). 


Hence, there are four cases as to where R 2 can be caught by R\ that found the exit. As before, we 
omit the extra cost 1 which is the time needed for the deployment phase during the case distinction. 


Case 1 : R 2 is caught in its pre-detour phase: The cost of Algorithm C(x, 0, A) is at most x + /(x)> 
exactly as in the analogous case of Algorithm £>(x, (p). 

Case 2: i ?2 is caught in its detour subphase-1: As in case 2 of the analysis of Algorithm £>(x, <P), if 
E is the position of the exit then, y = EA satisfies y >x~ /(x) (the relevant figure in this case 
remains the left-hand side of Figure 4). As longs as q(y) remains less than A, the cost of the 
algorithm remains y + 2 h(y). In order to find the maximum y for which this formula remains 
valid, we need to solve the equation A = q(y). This is possible by recalling that h(y ) = X+Q(y)~ 2b 
and by invoking the formula of h(y) as it appears in Lemma 4. By the monotonicity of h(y), 
we have that there exists unique ip satisfying h(ip) = X + A — ip. It follows that the cost of the 
algorithm in this case is at most sup x -f(x) <y< ^{y 
Case 3: R 2 is caught in its detour subphase-2: In this case, the relevant figure is the right-hand 
side of Figure 5. Let the exit be at point E, and let K denote the meeting point of the robots 
on the line segment GC. We set h'(y) := EK, which is calculated next in Lemma 5 (a). We 
conclude that in this case the cost of the algorithm is at most sup ^< y < x {y + 2h'(y)}. 

Case 4: R 2 is caught in its post-detour phase: Let td again be the total time a robot needs till it 
enters its post-detour phase. As in case 3 of Algorithm £>(x, </>), the cost of Algorithm C(x0, A) for 
this case is at most td+2p(x)- It thus remains to show how to calculate td — AB +BG+GC+CB, 
which is done in Lemma 5 (b). 

Lemma 5. The following statements hold for Algorithm C(x, 0, A): 

(a) Suppose that R\ finds the exit and meets R 2 in its detour subphase-2. Then the time h'(y ) that 

R\ needs from finding the exit until meeting R 2 is h'(y ) = 7V(x, 2 /, A, cp)/D(x 5 V, A, (p), where 


N(x, yAA)-= 2 + A 2 + (A + X - y) 2 + 2A(sin (</> - x) - sin (0 + y)) 

+ 2(A + x - y)(sin(x) + sin (y) - A cos (0)) - 2 cos (x + y ), 
D{x,y,\4>) :=2(A + X - y + sin (x) + sin (y) - A cos (</>)). 


(b) Suppose that R\ finds the exit and meets R 2 in its post-detour phase. Then, the total time that 
i ?2 spends in its detour phase is 











Proof. As an illustration of the proof we refer to Figure 6, which is a continuation of Figure 5. Let 
H be the symmetric point of E with respect to diameter AA'. As in Figure 4, L is the projection of 
H onto the supporting line of DB , and 6 denotes the angle BHL , whose value is given by (5). Also 
G r and C’ are the projections of G and C, respectively, onto DL. The calculations below follow the 
spirit of the arguments in Lemma 4. 




H 


L 


Fig. 6. The analysis of Algorithm C(x, 4>, A). 


(a) As before, y denotes the distance of the exit from point A. We have that EK — eA + hA + 
and therefore 


( 8 ) 


EK = EH + HB + BG + GC + 

2 (k 3 -H^ + E^-B^ + E^-Gd + H^-B^ + H^-Gd + B^- , 

where EK = h'(y), ~EH = 2 sin (y), ~HB = 2 sin (*=»), ~BG = A, and GK = q'{y). The inner 

products EF} ■ nh, Eli ■ BG, Hli ■ B<G are calculated exactly as in Lemma 4. For the remaining 
inner products we see that 

Eti ■ Git = ~EH GK cos (tt) = -~EH GK, 

-G^ = TTB GK cos = HBGK sin ((9), 

B& ■ Git = Wi GK cos (</>). 


Substituting the above in (8), we obtain an equation for h'(y) as a function of q f (y), y , X? an( i ( t ) - 
In the latter equation we can substitute q'(y) using the meeting condition (7), according to which 
q'(y) = y + h'(y) — x ~ A. Resolving the resulting equation for h'(y) gives the desired formula. 


(b) We need to calculate BG + GC + CB , where BG = A. First, we obs erve that GC = G'C' = 
BC' — BG' = sin (x) — A cos (</>). In order to calculate CB , we see that bA = bA + g3. Hence we 
obtain 

BC 2 = BG 2 + GC 2 + 2 ■ G& 

= A 2 + (sin (%) — A cos (</>)) 2 + 2A (sin (%) — A cos (0)) cos (0) 

= sin 2 (x) + A 2 sin 2 (0), 




































which concludes our claim. 


Before stating our main theorem, we summarize the total time required by Algorithm C(x, 0, A) 
in the following lemma. 


Lemma 6. The cost of Algorithm C(x, 0, A) can be expressed as 


(pre-detour phase )' 
(detour subphase-1) 
(detour subphase-2) 


' x + fix) 

s up x -f( x )< y <^{y + 2 h(y)} 

1 + max <( sup ip<y< x {y + 2h \y)} 

x + A + sin (x) - A cos (<f>) + 0in 2 (x) + A 2 sin 2 (</>) + 2 p(x) 

(post-detour phase) > 

where the functions f and p are as in (1) and (3), respectively; functions h(y) and h'(y ) are 
expressed explicitly in Lemmas 4 5 (a), respectively; and if is the unique solution to the equation 

MVO = x + a-y>. 

Using the statement of Lemma 6 and numerical optimization, we obtain the following improved 
upper bound. 

Theorem 1. For xo = 2.631865,00 = 0.44916 and Ao = 0.05762, the evacuation algorithm 
C(x o 5 0o?^o) has cost no more than 5.628. 

Proof. We examine the cost of our algorithm depending on where the meeting point of the two 
robots occurs. The guidelines of the analysis are suggested by Lemma 6. Also the deployment cost 
of 1 will be added at the end. Any calculations below are numerical, and were performed using 
MATHEMATICA. 

For the given parameters, we see that /(xo) — 1.99603, p(xo) — 0.506932, yo — f(x o) — 0.63584, 
and 0 = 0.755204. If the meeting point is during the pre-detour phase, then the cost is Xo + /(xo) < 
4.62791 (note that xo < x o)- If th e meeting point is in the post-detour phase, then the cost is 

Xo + A 0 + sin (xo) - A cos (<f> 0 ) + 0in 2 (xo) + Ag sin 2 (0 o ) + 2p(xo) < 4.627965. 

For the more interesting intermediate cases, we see that 

-0.5 y 2 + 3.45042?/ + a(y) + b(y) - 6.61768 


Kv) = 
h\y) = 


y - 0.900812 sin(y) - 0.434209 cos(y) - 3.45042 

—0.5 y 2 + 3.12552 y + (y - 3.12552) sin (y) - 0.847858 cos(y) 

y — sin(y) — 3.12552 


5.74387 


where a(y) (0.900812?/ — 2.85875) sin(y) and b(y) := (0.434209?/ — 2.01566) cos(y). We can see 
then that the cost in the detour subphase-2 is 

sup {y + 2 h(y)} < 4.627972, 

0.63584<y <0.755204 

and the cost in the detour subphase-3 is 

sup {y + 2 h'(y)} < 4.627961. 

0.755204<y <2.631865 


This completes the proof of the theorem. 








3 Lower Bound 


In this section we show that any evacuation algorithm for two robots in the face-to-face model 
takes time at least 3 + | + y/3 ~ 5.255. We first prove a result of independent interest about an 
evacuation problem on a hexagon. 

Theorem 2. Consider a hexagon of radius 1 with an exit placed at an unknown vertex. The worst 
case evacuation time for two robots starting at any two arbitrary vertices of the hexagon is at least 
2 + \/3. 

Proof. Assume an arbitrary deterministic algorithm V for the problem. V solves the problem for 
any input, i.e., any placement of the exit. We construct two inputs for V and show that for at least 
one of them, the required evacuation time is at least 2 + \/3. First, we let V run without placing an 
exit at any vertex, so as to find out in which order the robots are exploring all the vertices of the 
hexagon. We label the vertices of the hexagon according to this order (if two vertices are explored 
simultaneously then we just order them arbitrarily). Let t be the time when the fifth vertex, ^ 5 , of 


^X 




Fig. 7. Vertices of the hexagon as visited by algorithm V; t is the time when the fifth vertex V 5 is visited by some 
robot, say R±. One of the vertices adjacent to v$ has not been visited yet by a robot. 


the hexagon is visited by some robot, say i?i, i.e., robot R\ is at vertex v$ at time £, and four more 
vertices of the hexagon have been already visited, (see Figure 7). In other words, v$ and vq are 
the only vertices 7 that are guaranteed to not have been explored at time t — e, for any sufficiently 
small e > 0. Note that we must have t > 2, since at least one of the two robots must have visited 
at least three vertices by this time (and hence must have walked at least the two segments between 
the first and the second, and between the second and the third vertices visited in its trajectory). 

The first input I\ we construct has the exit placed at vertex vq. If t > 1 + \/3, then this 
input gives an evacuation time of at least 2 + y/3. Indeed, until robot R\ reaches P 5 , the algorithm 
V processes I\ identically to the case where there is no exit; further, at time £, robot R\ needs 
additional time at least 1 just to reach the exit. 

Hence assume that 2 <t< l + \/3- Let v x , v y , and v z be the three vertices that are non-adjacent 
to ^5 in the hexagon (see Figure 7). Note that the minimum distance between v$ and any of v x , 
v y , and v z is at least \/3. If vq £= {vx,Vy,Vz} then on input / 1 , V needs evacuation time at least 
t + y/3 > 2 + \/3, as R\ still has to reach the exit. 

7 It might be that V 4 and V 5 are explored simultaneously, or that V 5 and vq are explored simultaneously. In the 
former case ^6 is explored strictly after V 5 while in the latter V 4 is explored strictly before ^ 5 . 





Therefore, assume that ^6 0 { v x , v y , v z } and hence {v x , v y , v z } C {^ 1 , ^ 2 , ^ 3 , ^ 4 }. Note that, since 
t < 1 + >/3, on input p, robot Pi has visited at most one of v x , v y , and v z at time t. Hence, the 
other robot, P 2 has visited at least two of them. For the second input I 2 that we construct, we 
place the exit on the vertex v* that is the last vertex among v x , v y , and v z in the visiting order 
of the vertices by R 2 . Let t* be the time when R 2 reaches i?*, and note that at least until £*, the 
algorithm V behaves identical on the two inputs I\ and I 2 . As R 2 has visited at least one vertex 
before visiting v*, we have t* > 1 . 

Next we claim that R\ and R 2 cannot meet between time t* and t. The claim below shows that 
the former is impossible. Namely, we can prove: 

Claim 2 If t < 1 + \/3, then on input I 2 , Pi and R 2 do not meet between time t* (defined in the 
preceding paragraph) and time t. 

Proof. (Claim 2) Assume on the contrary, that on input p, Pi and R 2 do meet at some time t' at 
point P, with t* < t' < t. Observe that on input / 2 , robot Pi continues until time t' as on input R 
but having met P 2 at time t' might continue differently after time t'. Let ts = t' — t* be the time 
that P 2 uses on input R to get from the exit v* to P, and let tA — t — t' be the time that Pi uses 
on input R to get to vertex v§ from P. As v* and v$ are at distance at least \/3, and since t* > 1, 
we have 

y/3 < tA + ts = t — t* + t f — t* = t — t* < t — 1. 

So we obtain + 1 < t, which contradicts the assumption t < 1 + y/3. This proves the claim. ■ 

Having proved the claim, we conclude that on input / 2 , Pi continues until time t as on input 
R. Hence Pi needs at least t + \/3 > 2 + y/3 time to reach the exit on input R. This completes the 
proof of the theorem. ■ 

It is worth noting that the lower bound from Theorem 2 matches the upper bound of evacuating 
a regular hexagon, when the initial starting vertices may be chosen by the algorithm. Consider a 
hexagon ABCDEF and suppose that the trajectory of one robot, as long as no exit was found, is 
ABDC. Similarly, the other robot follows the symmetric trajectory FECD ; cf. left-hand side of 
Fig. 8. By symmetry it is sufficient to consider exits at vertices A , B or C. An exit at C is reached 
by each robot independently, while both robots proceed to an exit at A or B after meeting at 
point M, the intersection of segments BD and EC. Altogether, they need a total time of at most 
max{l + 4/^/(3), 1 + (2 + \fl)/\/ 3, 1 + y/3 + 1} to evacuate from the hexagon. It is easy then to 
verify that, in each case, the evacuation time of this algorithm is always upper bounded by 2 + \/3. 

In the above algorithm, the robots meet at M, regardless of whether the exit has been already 
found or not. The idea of our algorithm for disk evacuation presented in the previous section was 
influenced by this non-intuitive presence of a forced meeting. 

Combining Theorem 2 with some reasoning from measure theory, we obtain the following lower 
bound for our evacuation problem. 

Theorem 3. Assume you have a unit disk with an exit placed somewhere on the boundary. The 
worst case evacuation time for two robots starting at the centre in the face-to-face model is at least 
3 + f + 5.255. 



Fig. 8. The trajectories for R\ (red) and R 2 (blue) for the hexagon evacuation algorithm having evacuation time 
2 + >/3, while the exit has not been found, are depicted on the left. Right-hand side: At time 1 + — e, there is 

regular hexagon all of whose vertices are unexplored and lie on the boundary of the disk. 


Proof. It takes 1 time unit for the robots to reach the boundary of the hexagon. By time t — 1 + f, 
any algorithm could have explored at most ^ of the boundary of the disk. Hence for any e with 
0 < e < t, there exists a regular hexagon with all vertices on the boundary of the disk and all 
of whose vertices are unexplored at time t — e; see the right-hand side of Figure 8. Now, invoking 
Theorem 2 gives the bound of at least 1 + | + 2 + v^3 to evacuate both robots. ■ 


4 Conclusion 

In this paper we studied evacuating two robots from a disk, where the robots can collaborate using 
face-to-face communication. Unlike evacuation for two robots in the wireless communication model, 
for which the tight bound 1 + + \/3 is proved in [9] , the evacuation problem for two robots in the 

face-to-face model is much harder to solve. We gave a new non-trivial algorithm for the face-to-face 
communication model which improved the upper bound in [9]. We used a novel, non-intuitive idea 
of a forced meeting between the robots, regardless of whether the exit was found before the meeting. 
We also provided a different analysis that improved the lower bound in [9] . 

We believe that none of our bounds are close to be tight. More specifically, we do know that our 
upper bound is not optimal, since by disallowing robots to meet without having found the exit (by 
slightly truncating their trajectory), we can provably improve the performance of our algorithm. 
Unfortunately, the improvement we obtain this way is negligible (affecting the third significant 
decimal digit) while the additional required technicalities would be overwhelming, without offering 
new insights for the problem. This also suggests that the choice of the parameters we choose for 
our algorithm are not optimal. We are also certain that the proposed algorithm, i.e. family of 
trajectories we consider, cannot give the optimal trajectory, as it is intuitive that the optimal 
solution should be related to a properly defined differential equation ensuring that if robots meet 
during the deployment phase then the overall cost stays constant. Similarly for the lower bound, 
we believe that our proposed technique will serve as a guideline towards a more refined analysis 
that would reduce the gap. To conclude, a tight bound still remains elusive. 
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